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The theory of statistical mechanics is studied in the presence of Lorentz-violating background 
fields. The analysis is performed using the Standard-Model Extension (SME) together with a Jay- 
nesian formulation of statistical inference. Conventional laws of thermodynamics are obtained in 
the presence of a perturbed hamiltonian that contains the Lorentz violating terms. As an example, 
properties of the nonrelativistic ideal gas are calculated in detail. To lowest order in Lorentz viola- 
tion, the scalar thermodynamic variables are only corrected by a rotationally invariant combination 
of parameters that mimics a (frame dependent) effective mass. Spin couplings can induce a tem- 
perature independent polarization in the classical gas that is not present in the conventional case. 
Precision measurements in the residual expectation values of the magnetic moment of Fermi gases 
in the limit of high temperature may provide interesting limits on these parameters. 



I. INTRODUCTION 

The notion that the minimal standard model serves 
as a low energy limit to a more fundamental theory 
which includes a quantum description of gravity has led 
to the development of theories which extend the stan- 
dard model and predict qualitatively different physical 
phenomena Q. In this context, a framework for study- 
ing the effects of spontaneous Lorentz symmetry break- 
ing and possible resulting CPT violation within the con- 
text of conventional quantum field theory has been devel- 
oped 0, 13 and studied intensively . This framework 
formally contains all possible operators utilizing stan- 
dard model fields that satisfy coordinate reparametriza- 
tion invariance and is called the Standard-Model Ex- 
tension (SME). The minimal version of the SME is re- 
stricted to exhibit a number of useful properties: it pre- 
serves energy-momentum conservation, observer Lorentz 
invariance, conventional quantization, hcrmiticity, micro- 
causality, positivity of the energy, gauge invariance and 
power counting renormalizability. In addition, the frame- 
work provides a generic model for any theory which ex- 
tends the standard model and provides for spontaneous 
Lorentz symmetry breaking and CPT violation. For ex- 
ample, Lorentz violation may arise in noncommutative 
field theory p| or random dynamics models [||. 

To date, there has been no confirmed evidence for 
Lorentz violation, hence it is reasonable to assume that 
any violation must be small in conventional laboratory 

tperi: 
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have reached a precision that probes these parameters at 
Planck-suppressed scales. Recent analysis have also been 
extended to the neutrino sector pjj , |2J| , instantons [25j , 
supersymmetric models |26j| , and the gravitational sector 
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Statistical mechanics has been applied to rotationally 
invariant CPT-violating terms in the SME to provide a 
mechanism for baryogenesis in thermal equilibrium |28| . 
but a detailed exposition of the formalism for general 
Lorentz violation is lacking. It is the goal of this paper 
to address this issue. 

A general framework is provided in this paper for per- 
forming statistical mechanics calculations within the con- 
text of the SME. The required assumptions regarding 
statistical inference and the definitions of relevant ther- 
modynamic quantities are given. A complete analysis of 
the effects of all Lorentz-violating terms on a nonrela- 
tivistic ideal gas is performed. The corresponding ther- 
modynamics is studied and the various changes in ther- 
modynamic quantities are analyzed. In section II, the 
notation and the relevant formalism for performing sta- 
tistical mechanics calculations within the context of the 
SME is presented. Section III contains calculations of 
the relevant thermodynamic quantities for single particle 
systems. Section IV generalizes the result to a classical 
gas with a variable particle number and incorporates the 
chemical potential. Section V generalizes the result to 
a quantum gas of fcrmions and section VI gives corre- 
sponding results for spin-0 bosons. The conclusions are 
gathered in section VII. 



II. NOTATION AND FORMALISM 

Let {4>i}'iLi be a collection of states and suppose that 
{fj}j—i is a finite collection of real valued functions on 
the collection of states. Given a distribution of states, 
Qi = 'Z(V'j)' we denote by brackets the corresponding ex- 
pectations; (fj) = We denote the (informa- 
tion) entropy associated to the distribution qi by 



S = -fcy^gjlnfa) 



(1) 
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where A: is a positive constant that will later be identified 
with Boltzmann's constant. 

Suppose we are given an observation of the mean val- 
ues {(fj)}j—i and we seek a distribution that best fits 
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our observation. In his approach to statistical inference, 
Jaynes [2!| argues that one should choose a distribution 
which maximizes Q subject to the constraints on expec- 
tations given by the observations. In so doing, a now 
standard argument via variational calculus leads imme- 
diately to the solution 



z 



(2) 



where the Lagrange multipliers, Aj, are real con- 
stants and Z is the partition function, Z(X) = 

^. e -Ei =1 A '/3(^). prom © it follows that 
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If, in addition to being functions of state, the fj are per- 
mitted to depend on a finite collection of parameters, 
a ji 1 < j < m, then it is easy to see that 



d 
da k 
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The resulting maximal entropy is given by 
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(4) 



(5) 



Variation of the maximal entropy with respect to the 
parameters Xj, aj, and (fj) yields the relation 



(6) 



As is clear from our development, are formal 

rules of inference which follow for any system described 
as above, under the assumption of maximum entropy. 
Jaynes applies these formal rules to the study of equi- 
librium thermodynamics. More precisely, suppose that 
we consider a system of identical particles constrained to 
lie in a box of volume V. Suppose that f\=E denotes 
the energy levels corresponding to the possible states of 
the system and that fi = N denotes the number of par- 
ticles in the system. Given observations of the mean 
values, (E), (N), one obtains the partition function for 
the grand canonical ensemble: 



Z(a,P)=J2e-^- aNi 

i,3 



(7) 



where we have identified the Lagrange multiplier Ai with 
the scaled inverse temperature (3 — and the Lagrange 
multiplier A2 = a with — /3/i, where \x is the chemical 
potential. 

Allowing the energy to depend on the volume V and 
setting a\ — V in (@J we get the usual description of 
pressure 



P = 



dE\ 18 



(8) 



Equations GJ-© gi ye the central features of thermody- 
namics including expressions for the entropy as a func- 
tion of temperature, volume and particle number, and 
the usual first law of thermodynamics. The standard ex- 
pression for average energy and particle number are given 
by (JjJJ and the expressions for specific heat at constant 
volume and at constant pressure follow by taking the ap- 
propriate derivatives of the entropy: 



CV = T 



dS_ 
df 
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(9) 



(10) 



Thermodynamic potentials such as the Hclmholtz free 
energy are defined in the usual way. 

As was emphasized by Jaynes, the method of maxi- 
mum entropy provides accurate predictions of thermo- 
dynamic properties, assuming empirically accurate mean 
values and laws of motion embodied in the associated 
hamiltonian. For the purpose of probing the thermody- 
namics of systems of particles with symmetry violations, 
the method offers a simple framework for generating re- 
sults. All that remains is to specify the relevant under- 
lying hamiltonian. 

In the case of free spin-i Dirac fermions jb of mass m, 
the SME is determined by the lagrangian |2| 



(11) 



where 
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In the above expressions a M , 6 M , c^, d^, e M , / M , gx^ 
and H^v are real fixed background parameters which 
determine the Lorentz violation. For low energy ap- 
plications the associated nonrelativistic hamiltonian H 
has been obtained using a generalized Foldy-Wouthuysen 
transformation To second order in p/m, this hamil- 
tonian is given bv|3l| 



with 



H = ^ + H LV 
Am 



llix = .1 • n,o' ■ C,'— ■ D jk ^-a k 
m m 

j? PjPk , r PjPkl 

*jk-^ H (jjki— — a 

2m 



(12) 



2m 



(13) 



where common terms in the original lagrangian l|ll|) are 
collected: 



A = (ao - tocqo - me ) 



(14) 



B, 



(-bj + mdjo - \me jk ig k i + \e jk iH k i) , (15) 
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Cj = [aj - m(coj + Cj ) 



(16) 



Djk = [—boSjk + m(d kj + d 00 Sjk) 

+meki m (^g m ij + gmoofiji) + £jklHlo\ , (17) 



Fjk = -2 [(cjk + ^c 00 5 jk )] 



(18) 



Gjki = 2{[(d, 



- \{b.j/m + d jQ + ^e jmn (g mn0 + H mn /m))] 5 U 
+ \{h/m+ ^ei mn g mn0 )Sjk 

— ejlm(gmOk + Srofco)} • (19) 

The hamiltonian is invariant under spatial translations, 
therefore the logarithm of the associated grand partition 
function scales with volume as in the conventional case. 
Using Eq. (JSJ together with this property yields the re- 
lation 



PV 
~~kT 



HZ) 



(20) 



III. SINGLE PARTICLE SYSTEMS 



We first consider a system consisting of a single free 
spin-i particle governed by the hamiltonian H appear- 
ing in (|12p. constrained to a cube of side length L. For 
small violations of Lorentz symmetry, we derive the cor- 
responding perturbations to the statistical mechanics us- 
ing the single particle partition function. This simplifies 
the initial analysis by eliminating the need to discuss the 
chemical potential. 

The standard unperturbed solutions are written in the 
form 



<iw = ri sin 



rii'KXi 
L 



(21) 



where n = (ni, n^) is a triple of positive integers and 
s G {1,-1} denotes a sign. Note that the two-component 
spinor Xs may depend on n. The corresponding unper- 
turbed energy levels are written as 



n ' s 2mL 2 



(22) 



The first order correction to the energy levels due to 
the Lorentz- violating terms are found using standard de- 
generate perturbation theory as: 



TL.Sf 



2mi 2 



{An 2 + ^F lin j + s\g(n)\), 

(23) 



where the vector £7(n) is defined with components 
2mL 2 



ir 2 ?i 



(24) 



The n dependent spinors Xs have been chosen to satisfy 
G(n) ■ u\ s = s Xs an d therefore diagonalize H LV . Note 
that C and D perturbations depend linearly on momen- 
tum and therefore don't contribute a correction to the 
energies. 

The perturbed energy is written as 



EL 



E$?l + 6E„ 



(25) 



with 5E n ^ s given by the matrix elements (|23[1 . The par- 
tition function for a single particle becomes 



-0E„ 



(26) 



Approximating the sum on the right hand side of i|26|) 
by the appropriate integrals, we obtain expressions for 
partition functions corresponding to hamiltonians with 
Lorentz-violating terms. The partition function is even in 
the spin-dependent correction terms and therefore there 
are no lowest order corrections contributed by B or G. 
The A and F terms correct the partition function as 



2e-f> A n Q V(l - ±Tr(F)) , 



(27) 



where V is the volume of the box, uq = (m/2/3TTh 2 ) 3 ^ 2 
is the quantum concentration, and Tr(F) = J^Fu- 

Using the relationship (J5W) = -(d/dj3)ln{ZW) (cf 
©), it follows from (|27[) that only the A term corrects 
the energy. The A term corresponds to a constant shift 
in all of the energy levels, hence it is possible to rede- 
fine the zero point of the energy |32l | to eliminate this 
contribution. 

The correction to the partition function due to the F 
term can be incorporated into an effective mass for the 
fermion 



= (1 - ±Tr(F))r 



(28) 



because only the rotationally invariant trace component 
contributes. The Tr(F) term appears to be a trivial scal- 
ing that can be absorbed into the mass, but this is not 
quite true because the effective mass violates boost in- 
variancc. This means that the effective mass can be dif- 
ferent in various laboratory frames. However, such effects 
will be relativistically suppressed and therefore difficult 
to observe using Earth-based experiments. 

More interesting is the correction to the expectation 
value of the spin. In this case, only the B and G terms 
will contribute. The expectation is calculated as 



(s (1) ) = (^rEjf^-^- 2 ^ 

2mL 



(n)|) 



^- X /^£0(n)e-^ 



-/m--Tr(G) , 



(29) 
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where the vector (Tr(G))fc = J2iGuk is defined. Note 
that while it is not in general possible to simultaneously 
diagonalize the contribution from B and G terms, the 
same computation allows us to treat the case where both 
terms occur. The corresponding spin expectations de- 
couple. 



IV. CLASSICAL GAS 

We now introduce the chemical potential and consider 
a classical gas of free spin-i Dirac Fermions. As in sec- 
tion II, N denotes the (variable) particle number for the 
system. The grand partition function for the system can 
be written in terms of the single-particle partition func- 
tion as 



(30) 



where — —ft^ ' and fi^ c ' is the chemical potential 
of the classical gas. Thus equation l|27l) also gives the first 
order corrections for the grand partition function. Using 
the grand partition function, the resulting expressions for 
expected particle number and energy are 



<£(C)) = t( N (c) )kT 



(31) 
(32) 



bmce PV = kT\nZ {c '\ it follows that there is no change 
in the classical ideal gas law, even in the presence of 
Lorentz violation. 

Solving for the chemical potential ijS c ' (to lowest order 
in Lorentz violating parameters) yields 



t ( c > = -kT In 



2ng 

n(C) 



- |Tr(F) 



(33) 



where nq is the quantum concentration and n^ c ' = 
(N^)/V is the concentration of the classical gas. We 
conclude that in the presence of Lorentz violation there 
is a change in the chemical potential as expected from 
the effective mass argument presented in the previous 
section. 

The corresponding entropy is found using equa- 
tion JSJ) as 



= a ^(N^)k + ^(N^)k . (34) 
Using H33f) . the modified Sackur- Tetrode equation is 



= {N (c) )k 



(35) 



From (|35|l and (j9"|)- (|10[l it is clear that there is no change 
to the specific heat. 

Finally, we verify that the expectation of the spin is in 
fact the single particle result in equation (|29ll times the 
expected particle number 



(s( c >) = -(JV( c >)[/3B + iTr(G)] 



(36) 



V. QUANTUM GAS - FERMIONS 

Next, the quantum occupancy of the orbitals is re- 
stricted to be or 1 and the low-temperature limit is 
analyzed. With notation from previous sections, the par- 
tition function for the grand canonical ensemble associ- 
ated to a Fermi gas is 



Z«>(a) = J!! 1 



(37) 



In the unperturbed case, the grand partition function is 
calculated using the appropriate integral approximation 
as (zero subscripts represent unperturbed quantities) 



In 



(4 Q) K) 



n 2 In f l + er" ' ^ 



dn 
(38) 



where A = h/(2irmkT)2 is the thermal wavelength (re- 
lated to the quantum concentration by 1/A 3 = uq) and 
f v (e~ a ) is the Fermi-Dirac integral 



1 



T(v) J e a e x + 1 
Using J2Jl gives the standard results 

<iV (Q) (a )) = Jr^/f(e- ao ) , 



dx 



(39) 



(40) 



<4 Q) (ao)> = ^ iV o Q) («o))fcT^z^ , (41) 



and the quantum ideal gas law 

P fi(e~ a °) 



ikT 



h(e-*o) 



(42) 



As in the classical case, there are no first order cor- 
rections to the partition function for Lorentz violating 
terms given by the coefficients B, C, D and G appear- 
ing in i|13|) . For violating terms of type F, a change of 
variables in the integral (|38|l gives first order corrections 
for the partition function which we write as 

ln(z«>(a)) ~ (l-±Tr(F))ln(z^(a)) .(43) 

As in the classical case, only the rotationally invariant 
component of F contributes, therefore it is possible to 
absorb the term into an effective mass as before. Stan- 
dard calculations then immediately give the results 

= (l-iTr(F)) A /|(e -«) , (44) 
= |(JvW)(a))fcT^i^ , (45) 
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as well as the ideal gas law 



P 



niQ^kT fs.(e- a ) 



(46) 



where n«) is the concentration of the quantum gas (not 
to be confused with the quantum concentration hq de- 
fined earlier). Note that the ideal gas law is in fact 
modified in the quantum case due to its dependence on 
q ^ oq. The correction due to F may be incorporated 
easily by replacing A(m) by A(m*) using the effective 
mass given in equation l(2"81) . 

From (|39|l it is clear that the map a — ► /|(e _a ) is 
strictly monotonic and thus invertible. Writing the in- 
verse as T, we have a formal expression for the chemical 
potential 



(Q) 



-kTT 



(47) 



Equation 1)47(1 demonstrates that conventional formu- 
las can be used to obtain the relevant quantities. These 
formulas can be found in a standard statistical mechanics 
text such as [s^. For instance, the chemical potential at 
zero temperature defines the Fermi energy and is modi- 
fied as 



M «)(T = 0) = ^~£^ o >(l + iTr(F)) , (48) 

where Ej^ 1 = (h 2 /2m)(3w 2 n ( -^) 2 / 3 is the conventional 
Fermi energy. Approximating to the next highest order 
in temperature gives 



12 



kTY 



(49) 



Using i|44p . 145|) . an asymptotic expansion for fs and 
(g5J) we obtain 



(N(Q)) 



5tt 2 /fcT 
"12" ^ 



(50) 



From © and 1(51) [) we get an expression for the specific 
heat in the limit of low temperature 



C v 



{N(Q))k 



it 2 kT 



7T 2 kT 

~2~^C°) 



(1 - iTr(F)) 



(51) 



Writing the entropy using J^J gives the perturbed entropy 
in the low temperature limit as 



(52) 



as expected from the integrated equation for specific 
heat. The low-temperature ideal gas law perturbation 



5 * 



5^ 
~L2~ 



jfeTV 



(53) 



The expectation value for the spin can be calculated 
in the quantum regime using the fractional occupancies 



f(n,s) 



1 



(54) 



where E n-S — E n . s + 5E n ^ s and the perturbation is as 
given in 123fl . With notation as above, the expectation 
of the spin is given by 



\S(n) 

7T 2 h 2 

2mL 2 

{N (Q) } 



-2(3 



2rn. 1 ,2 



(1 



o 



2^f 1 _(e- a )B + ^Tv(G) 



(55) 



This calculation demonstrates the surprising fact that the 
contribution of the G term is temperature independent 
and does not randomize at high temperature. At low 
temperatures, the contribution from the B term can be 
written as 



3 B 



i-^y 

12 \E F ) 



(56) 



This is identical to the result for a quantum Fermi gas 
in an external magnetic field when only spin interactions 
are relevant. Note that the parameter B is a fixed back- 
ground vector that does not rotate with the experiment. 



VI. QUANTUM GAS - BOSONS 

It is possible to generate a model for a free spin-0 boson 
gas by combining two fermions into a singlet representa- 
tion of the spin group. This means that the expectation 
value of all spin-couplings vanish in the hamiltonian 113(1 . 
The resulting hamiltonian is given by 



H 



P 

2m 



A + C. 



Pi 



■ F, 



PjPk_ 
2m 



(57) 



Choosing the ground state energy to be zero and employ- 
ing the notation of the previous sections, the associated 
grand partition function for the unperturbed case is 



Hzf B \a )) = - Vln(l-e- a °e-^) 



ln(l - e - Q0 ) 



(58) 



where E^ — (Tr 2 h 2 /2mL 2 )n 2 as before, and the ground 
state has been separated out to allow for Bose-Einstein 
condensation at low temperatures. Approximating the 
sum as an integral gives the standard result 

ln(Z« B )(a )) = ^^5|(e- Q0 )-ln(l- e - Q0 ) , (59) 
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where A is the thermal wavelength and ff|(e a ) is the 
Bose-Einstein integral 

9u{e- a ) = — -r / ^~^dx . (60) 
r W Jo e e - 1 

Using @ we obtain the expected number of particles in 
the excited states and the associated energy as 

(iV (QS) (a )) - (Nao) = ^V 9l (e-^) , (61) 
(4 Q) («o)> = |^^3|(e- ao ) , (62) 

where (Ngo) = [e a ° — is the expected number of 
particles condensed into the ground state. 

The only nontrivial leading order perturbation in l|57|) 
arises from the F term. A calculation which follows that 
done for the case of fermions gives 

ln(z« fl )(cO) * (l-iTr(F))^U 5| (e- Q ) 

-ln(l-e- Q ) . (63) 

It follows from © and JU^-GHl tnat f° r the perturbed 
case we have 

(N^ B \a))-{N G0 ) = (l-i-ft(F))i ff| (e-«) , 

<£W B )(a)) = |fcT^ 5s ( c -) . (64) 

As in the Fermi case, the chemical potential can be ex- 
pressed as a function of the number of particles in ex- 
cited states. The computations, modulo obvious modi- 
fications, are similar to l|47|) - (|53l) . Because only Tr(F) 
enters into the grand partition function, it is possible to 
use the concept of effective mass to absorb the effect as 
before. Standard results of Bose-Einstein condensation 
therefore hold in a given laboratory frame. The effect is 
nontrivial, as in the fermion case, because boosting the 
experiment will change the effective mass. Modifications 
of these calculations are required if the ground state wave 
function does not exhibit rotational symmetry, such as is 
the case in actual condensate experiments due to some 
optical and magnetic trapping configurations [34[ . In this 
case, the zero point of energy can depend on orientation 
of the apparatus and induce rotational variations in the 
condensate properties. An analysis of this type for spe- 
cific experiments would be interesting, but is beyond the 
scope of the paper. 

VII. CONCLUSION 

The formalism of statistical mechanics in the pres- 
ence of symmetry violation parallels the conventional sit- 
uation. The laws of thermodynamics are the same as 
in the conventional case, although the specific expecta- 
tion values of thermodynamic quantities can be modi- 
fied by the Lorentz-violating terms. Our approach in- 
volves an assumption of a maximal lack of information 



(or entropy) regarding a system subject to various con- 
straints imposed by the physical observables. The tem- 
perature and chemical potentials can be defined simply as 
the lagrange multipliers associated with the constraints, 
hence two systems in equilibrium automatically have the 
same temperature and chemical potential as there is only 
one lagrange multiplier for each overall constraint. This 
method in fact produces equivalent results to the more 
conventional assumption of equal a-priori statistics, how- 
ever, it has the advantage of providing straightforward 
definitions of thermodynamic variables with conventional 
equilibrium properties. 

As an explicit example, the unperturbed system was 
assumed to be an ideal gas in the absence of any ex- 
ternal applied fields (such as magnetic or gravitational 
fields). Expectation values for scalar thermodynamic 
quantities such as energy and particle number were unal- 
tered except for an overall scaling factor Tr(F). This 
happens because the unperturbed system is rotation- 
ally invariant and scalar expectation values can only be 
corrected by perturbations with commensurate symme- 
try. This can also be incorporated as an effective mass 
to* = (1 — ^Ti(F))m in the hamiltonian, although the 
effective mass defined in this way depends on the ob- 
server's Lorentz frame. For instance, a gas in motion 
on the surface of the Earth would exhibit slight sidereal 
variations in effective mass due to changes in Tt(F). This 
effect is relativistically suppressed and unlikely to be ob- 
servable for physically reasonable values for the violation 
parameters. 

More interesting is the net spin expectation value con- 
tributed by the terms that couple to the spin. The pure- 
spin coupling Bj mimics a constant background magnetic 
field and induces a corresponding magnetic moment per 
unit volume in the gas. Any additional applied magnetic 
field could be added to this term to calculate the net 
result on an actual experiment. For example, an Earth- 
based experiment will rotate in space and the constant 
background vector will interfere with the applied mag- 
netic field to produce a sidereal variation. 

The derivative-spin coupling Gijk generates a funda- 
mentally new type of effect that induces a temperature- 
independent polarization in the classical gas that is pro- 
portional to Tr(G). This means that even at very high 
temperatures there will be a net residual polarization 
that does not randomize. This is a reasonable result, 
considering that the Tr(G) term couples the spin to the 
conventional kinetic energy term in the hamiltonian and 
indicates that any physical effects should scale accord- 
ingly. An effect of this type should be clearly distin- 
guishable from stray magnetic field effects that have an 
inverse temperature dependence. 

In the zero temperature limit, both the Bj and Gijk 
terms contribute a polarization to the Fermi gas. The B 
contribution depends on n 1 ^ 3 while the Tr(G) depends 
linearly on particle density n. This means that the effects 
of the Tr(G) term can grow significantly faster with in- 
creasing density than conventional magnetic field effects. 
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